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ABSTRACT 

We present a newly developed approach to solar coronal magnetic field ex- 
trapolation from vector magnetograms, based on the Principle of Minimum Dis- 
sipation Rate (MDR). The MDR system was derived from a variational problem 
that is more suitable for an open and externally driven system, like the solar 
corona. The resulting magnetic field equation is more general than force-free. Its 
solution can be expressed as the superposition of two linear (constant-a) force- 
free fields (LFFFs) with distinct a parameters, and one potential field. Thus the 
original extrapolation problem is decomposed into three LFFF extrapolations, 
utilizing boundary data. The full MDR-based approach requires two layers of 
vector magnetograph measurements on solar surface, while a slightly modified 
practical approach only requires one. We test both approaches against 3D MHD 
simulation data in a finite volume. Both yield quantitatively good results. The 
errors in the magnetic energy estimate are within a few percents. In particular, 
the main features of relatively strong perpendicular current density structures, 
representative of the non-force freeness of the solution, are well recovered. 

Subject headings: MHD — methods: data analysis — Sun: corona — Sun: magnetic 
fields 
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Introduction 



Plasma relaxation processes are ubiquitous in astrophysical and laboratory plasmas. 
One well-known approach is to invoke the variational principle, minimizing the t otal 



magn etic field energy subject to the constraint of constant total magnetic helicity (iTayloi 



19741 ) ■ Such an approach yields the magnetic field equation for a relaxed plasma state, in 
which the Lorentz force vanishes, i.e., the so-called linear force-free field (LFFF) with a 
constant a parameter, 

VxB = aB. (1) 

Although the above equation ([1]) has been widely applied in modeling magnetic field 
configurations in laborato ry plasma devices and in solar a tmosphere, it is found inadequate 



i n cer tain situations (e.g.. 



Gary 



2001 



Metcalf et al 



19951 ). In particular. 



Amari &: Luciani 



(I2OOOI ). among others, showed by 3D numerical MHD simulation that in certain solar physics 
situation, after initial helicity drive, the final "relaxed state is far from the constant-a linear 
force-free field that would be predicted by Taylor's conjecture" with helicity conserved. 
They suggested to derive alternative variational problem. 

An alternative approach, also based on the variational principle, was recently developed 
for an open system with external drive, and applied to theoretical investigation of the 



solar arcade structures with fiow f,Bhattacharyya fc Janaki 



2004 



Bhattacharyya et al 



20071 ) . It was later applied to develop a new approach to the extrapolation of solar 



coronal magnetic field in non-force free state (iHu fc Dasgupta 



2006 



2007 



Hu et al. 



20071 ). The core of the approach is to construct a variational problem based on the 
Principle of Minimum Dissipation Rate (MDR), which states "The steady state of an 
i rreversible proc ess is characterized by a minimum value of the rate of entropy production" 



(jPrigogine 



19471 ). In most cases, the entropy production rate is equivalent to the energy 



dissipation rate. The basis is the set of generalized momentum balance equations, and 
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the system, with flow, is alway s in broadly-defined force-bala nced dynamic equilibrium, 
satisfying the MHD equations 



( iMontgomery Sz Phillips 



1988 



Bhattacharwa fc Janaki 



Dasgupta et al 



1998 



20041 ). Prior works based on MDR 



2002) had successfully shown that it 



was able to yield a pressure-balanced configuration, supporting a finite plasma pressure 
gradient found in plasma confineme nt devices. Recently a proof of MDR was rigorously 



sought by 3D numerical simulations fiShaikh et al. 



20071 ). We refer the readers to the above 



referenced literatures for detailed description and justification of the theoretical basis of 
MDR. 

In short, analogous to the principle of minimum energy, but for a more complex, open 
system like the solar corona, the MDR yields the following equati ons for the magnetic field 



B and fiow vorticity a; = V x v, y being the plasma fiq w velocity (iBhattacharyya fc Janaki 



2004 



Bhattacharyya et al. 



2007 



Hu fc Dasgupta 



20071): 



V X V X B + aiV X B -h 6iB = V^, 

V X V X cJ + X uj + b2UJ = Vx- 



(2) 
(3) 



Here, coefficients a and b are constants, and involve the parameters of the system. The 
right-hand sides of both equations are arbitrary, undetermined functions. 

In what follows, we first briefiy provide a heuristic justification of the MDR in Section [2l 
based on a simple but standard analysis. We then present a brief but comprehensive 
description of the full MDR-based coronal magnetic field extrapolation approach in 
Section [3l including a new test case study utilizing 3D MHD simulation data of a bright 
point region, highlighting the non-force free features of the solution. In Section HI we 
develop a practical approach that requires only one single-layer vector magnetogram and 
show that it achieves the same satisfactory results as the full approach for the same test 
case. In the last Section [5l we conclude and discuss the significance and limitations of our 
approach. 
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Heuristic Justification of the MDR 



In a generally resistive plasma, the Lundquist number, S, is defined as the ratio of 
the resistive time scale over the Alfven time scale, S = ^ (e.g. 



Ortolani fc Schnack 



19931). 



The Lundquist number scales with Lq/tjo, the ratio of the characteristic length over the 
resistivity of the system. Therefore the number S can be very large when Lq is large and/or 
plasma is highly conductive {rjo is small). For example, it can reach ~ 10^ — 10* in a hot 
fusion plasma, while for solar coronal structures it often exceeds 10^^. 



From standard procedures (e.g. 



Ortolani fc Schnack 



19931 ). the decay rates of global 



magnetic helicity, K, and magnetic energy, W, are given as follows (with the magnetic field 
being Fourier decomposed as B = ^,^bkexp(2k ■ r)): 



R 



dK 


27] 


~(U ~ 




dW 




dt 


s. 



V 
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Taylor envisioned the r elaxation process oc c urring as a result of small scale turbulence 



with S ^ 1. Following 



Ortolani &: Schnackl (jl993l ). at scale lengths for which the Fourier 



decomposition wavenumber is ~ S*^/^, the energy decay dW/dt is ~ 0(1). But for such 
scale length, the helicity decay rate dK/dt is ~ 0{S^^^'^) <^ 1. Thus, from this heuristic 
argument, we may expect that small scale turbulence dissipate energy at a greater rate than 
helicity. In the same way, we can show that dR/dt^ Sfc ^^^fc §°^^ ^ much faster 
rate that dK/ d t. Furthermore, a r ecent 3D time-dependent simulation of MHD fiuids proves 



this argument (jShaikh et al. 



20071 ). It is shown that the energy dissipa tion rate, / ■n j '^dV, 



decays the fastest towards a minimum state. We again refer readers to 
for a complete description of the 3D numerical simulation results. 



Shaikh et al 



fl2007h 
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3. The Full MDR-Based Approach to Coronal Magnetic Field Extrapolation 



The full MDR-based approach to coronal r aagnetic field ex t rapo 



Hu fc Daseuptal (120071 ): 



ation has been 



HuetaL 



fl2007h 



developed and tested against analytic models in 
based on the general equation yielded from MDR. The following equation results by taking 
an extra curl on both sides of Eq. ([2]) : 



VxVxVxB + aiV x V x B + 6iV x B = 0. 



(4) 



The aim is to solve the above equation in a finite volume utilizing boundary data. The key 
is that one exact solution to Eq. (jl]) exists and can be expressed as the superposition of three 
LFFFs. Each satisfies Eq. ([1]) with distinct a parameters. For the sake of completeness, we 
briefiy reiterate the procedures of the full approach in the following subsection. 



3.1. Procedures 



The exact solution to Eq. (jlj) is written 



B = Bi + B2 + B3 



(5) 



with V X B,; = ajBj, i = 1,2,3. Subsequently, one obtains (taking curls on both sides of 
Eq. ©) 

/ B \ /b,\ 



V X B 
\ V X V X B / 



V 



B2 
VB3/ 



where the matrix V is a (transposed) Vandermonde matrix composed of constant elements 
ai~^,i,j = 1,2,3. It is guaranteed inversible provided that the parameters ctj, i=l,2,3, are 
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distinct. Therefore, it follows 



B2 
VB3/ 



\ 



B 

V X B 
\ V X V X B J 



(6) 



The above equations provide the bottom boundary condition for each LFFF Bj, given vector 
magnetograph measurements B at certain bottom levels on the solar surface. In order to 
evaluate the double-curl term at the bottom boundary, ideally at least two levels of vector 
magnetograms have to be utilized and onl y the normal compone nt of Bj perpendicular 



to the bottom boundary can be obtained (IHu fc Dasgupta 



20071). Then each LFFF Bj 



is extrapol ated into a fin i te vol ume above the bottom boundary using the FFT-based 
method by 



Alissandrakid (jl98ll ) with known aj parameters for a Cartesian geometry. The 
summation of the three LFFFs yields the final solution according to Eq. ([5]), a non-force 
free field in general. 

It turns out that one of the a, has to be zero and the corresponding Bj becomes 



potential. The reaso r i is simp 



Hu fc Dasgupta 



2007 



Huetal 



y tha t the solution ([5]) has to satisfy Eq. 



see 



20071 . and Section H]). We arbitrarily choose 0:2 = 0, then 



— («! + 03) and hi = a^ai. 

The determination of ai and 03 is by tr ial-and-erro r: enumerate all possible pairs of ai 



and ^3, preferably within the limit of amax (IGary 



19891 ). for each pair and ^2 = 0, solve for 



the transverse field bj^ at the bottom boundary only, t hen evaluate the dev iation between 
the "exact" (measured) Bt and ht = J2^=i^it by (e.g., 

M M 



Schrijver et al. 



20061) 



En 



J]|B,,-b,,|/^|B 



t,i\y 



(7) 



i=l 



1=1 



where M = N"^, the total number of grids on the transverse plane, and the optimal pair 
(Q;i,a3)opt is chosen such that the corresponding E°p^ is minimum. 



- 8 - 



To summarize, the procedures of the full MDR-based approach to coronal magnetic 
field extrapolation constitute the following two steps: 

1. Set 0^2 = and search through LFFF parameters, ai and {ai > a^, due to 
the exchangeability of subscripts, 1 and 3). Apply an LFFF solver to calculate 
the transverse components bj^, utilizing the normal component Bin (given by 
equations ([6])), i = 1,2, 3, at bottom boundary only. A pair of optimal a^, i = 1,3, is 
found, for which the En, given by Eq. ([7]), is minimum. A plot of En distribution over 
the («!, as) parameter space is prepared to show the goodness and uniqueness of the 
solution. 

2. Solve for Bi and B3 in a finite volume above the bottom boundary, for the optimal 
«! and 03 found in Step 1, respectively, and the potential field B2 as well. Obtain 
B = Bi + B2 + B3 (Eq. ©). 

The run time of the approach based on the FFT algorithm solving each LFFF is 
approximately in the order of P ■ N"^ log2 A^, largely dependent on the number of search 
grids, P, on {ai,a^) in Step 1. However, this process can be easily parallelized by simply 
dividing the search domain into individual pieces and assigning each piece to a separate 
processor. One can also start with a coarse grid, then allocate a small but finer grid around 
a minimum, thus the total number of iterations can be effectively reduced. 



3.2. A Test Case Study Using Numerical Simulation Data 



Several analytic solutions have been utilized to test the full approach outlined in the 
previous subsection and the results showed the approach was able to recov er the solution in 



a finite volume to ce rtain degree of satisfactory accuracy for several cases (IHu fc Dasgupta 



2007 



HuetaL 



20071 ). Here in order to simulate of real magnetograph measurements 
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as closely as possible, we 
X-ray Bright Po i nt (BP) 



employ 3D MHD simula tion data surrounding a 



Otto et al 



20071 : 



Biichner 



Brown et al. 



2005 



200l|) from 



Biichner fc Nikutowskil (120051 ) (see also, 



regio n of an 



20061 ) ■ The simulation started with actual solar magnetic 
field measurements, and was continuously driven by photospheric motion inferred from 
solar observations. 

The data were provided in a finite Cartesian volume, x x y x 2;=128 x 128 x 63, with 
z being the vertical dimension along the normal direction, at one chosen moment during 
the 3D dynamic MHD simulation. The bottom two layers of vector magnetic field data are 
utilized to generate the bottom boundary [z = 0) conditions by Eqs. The two steps 
are carried out as described in Section [3m The distribution of En in (ai,a3) parameter 
space is shown in Figure [H The minimum value is E°^^ ^ 0.30, with corresponding 
{ai,a3)opt ~(0. 00156, 0.000779) (dimensionless), as denoted by the + sign in Figure [H 
Fairly significant amount of uncertainty associated with this optimal pair exists. However, 
further analysis shows that choosing any other pair of (ai, a^) within the innermost contour 
does not significantly chang e the extrapolati o n resu lts. Such similar pattern and small a 



values were also observed in 



Hu fc Dasguptal (120071 ). 



The corresponding optimal transverse magnetic field vectors ht at z = 0, are 
shown in Figure [2], together with the "exact" solution Bt from the simulation data 
and the corresponding potentia l field results using the normal component only 



( iVenkatakrishnan &: Gary 



19891 ). Our solution agrees very well with the "exact" one, 



whereas the potential field result exhibits apparent deviations. 

The final extrapolation results are shown in Figure [3l in the form of 3D field line plots. 
The bottom image shows the normal component distribution at ^ = 0, with two major 
polarities in the center corresponding to the location of BP. The field lines from our result 
demonstrate a good deal of similarity (nearly identical in the upper left part) to the "exact" 
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solution, indicating good-quality recovery of the magnetic field in the finite volume from the 



bottom two layers of data. This judgem ent is supportec 



so-called figures of merit, adopted from 



Schrijver et al 



jy th e set of quantitative measures. 



(120061 ) , given in Table [H They are 



used to evaluate, quantitatively, the agreement between the exact vector field B and the 
corresponding extrapolated field b in a volume (columns 2-5). The 2nd row in Table [1] 
lists the numbers for this calculation (result 1), and the last two rows list the values for a 
result identical to the exact solution and a potential field extrapolation, respectively. The 
last two columns show the energy estimates, the magnetic energy ratio of the extrapolated 
field b over the exact field B, i.e., e = b^/B^, and the energy ratio, e^, that of b over 
the corresponding potential field. For this case, several figures from the potential field 
extrapolation are comparable with our results. However, the energy estimates clearly mark 
the distinction: ours only differs from the exact ones by ~ —1%, whereas the potential field 
results differ by as large as ~ —16%. 

To further characterize the non-force freeness of the solution, which cannot be captured 
by either the potential or the force-free extrapolation, and to better visualize the results. 



we ca. 



[e.g. 



culate the fiel d-line integrated current density, J = f jdl, along individual field lines 



Biichner 



20061 ). Figure m shows such quantity separated into field- aligned (|| B) 



component, Jy, and the component perpendicular to B, J±, for the "exact" solution. A 
non-vanishing indicates the field is non-force free. Two footpoints of one single field 
line, both rooted on the bottom boundary, have the same value. Figure [5] shows the 
corresponding J\\ and J± for our result 1 discussed above. The distribution of J\\ is well 
recovered with concentrations at the two major polarities around the center. The major 
features of J± at the same locations are well recovered as well, but most of the remaining 
weak structures are lacking, compared with the right panel in Figure HI However, both J± 
and J|| enhancements near the lower boundary tracing back to one of the strong polarities 
to the right, remain. It is indicative of strong currents, implying the possible locations of 



current sheets that are important for magnetic reconnection (iBiichnerl 120051 . 120061 ). 



At presen t time, multip 



(but see, e.g., 



Metcalf et al. 



e (> 2) layers of vector magnetograms are hardly available 



19951 ). Only the photospheric vector magnetic fields are 



routinely observed. Sometimes, the chromospheric ma gnetic field line-of-sight (LOS) 



components were inferred (e.g.. 



Choudhary et al. 



20011 ). When the observations are near 



the disk center, the LOS component can be used as the normal component B^, then the 
term (V x V x B)^=-V^B^ may be approximated by (denoting the distance between the 
chromosphere and the photosphere by A^;), 

V^B, ^ VlB, + 2(5f ™™ - Bf"'" - —^Az)/Az^. (8) 

oz 

The first-order derivative dBz/dz is approximated by the divergence free condition using 
transverse magnetic field vectors at the photosphere. The 3rd row of Table [1] (result 2) lists 
the figures of merit of the extrapolation result using the bottom boundary magnetic field 
vectors and only the B^ component at one level immediately above to approximate V^-B^ 
by Eq. ([8]). It shows somewhat deteriorated results, and the error in energy estimate is 
~ +12%. 



4. A Practical Approach Utilizing Single-Layer Vector Magnetogram 

As mentioned earlier, the multi-layer vector magentograms are not routinely available. 
On the other hand, the high-quality photospheric vector magnetograms become increasingly 
available. In order to apply our method to actual measurements currently available, we 
adapt our approach to employ only one single layer vector magnetogram, at the expense of 
less-general, expected limited applications, as to be described below. 

From Eq. ([2]), considering the arbitrary nature of V'?/', one may write one exact solution 
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to ©, 

B = Bi + B3 + cBpot, (9) 

where a constant, c, and the potential field, Bpot, obtained from the known normal field 
at the bottom boundary, are introduced. Then it follows Vtp = bicBpot- This is equivalent 
to the full approach in Section [3l except that the potential field B2 is chosen a special form, 
cBpot, and is largely known subject to a undetermined constant. 

With an additional unknown parameter, c, the procedures are similar to those outlined 
in Section ISTTl but for a reduced, 2nd order system. The following equations, in place of ([6]), 
now provide the bottom boundary conditions for each LFFF from one single layer vector 
magnetogram: 

(a3-ai)5i. = a35;-(VxB)„ 
(ai - 03)532 = aiK - (V X B)^, 

with B'^ = Bz — cBpot,z- In addition to the (ai,a3) search domain in Step 1, another 
dimension in c is added. A desirable choice of allowable c values is c £ [—1, 1], to l imit t he 
dominance of the potential field. The approach first reported in IHu fc Dasguptal ( 120061 ) 
represents the special case, c = 0. 

We again test this practical approach against the same data set utilized in Section [31 
but only use the vector magnetic field data on the bottom boundary, z = 0. Only a limited 
number of c values are chosen, c G [—1, —0.5, 0, 0.5, 1]. The En distribution at Cmin = —1.0, 
which yields the minimum, E°^^ ^ 0.22, is shown in Figure [6l Compared with Figure [Tj 
the uniqueness of (ai,a3)opt is much improved. The corresponding at 2; = 0, and the 3D 
field line plot are shown in Figures [7] and [H respectively. The corresponding figures of merit 
are given in the 4th row (result 3) of Table [H Good agreement with the exact solution 
is achieved. The energy estimate has a ~ —5% error. The field-line-integrated current 
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densities shown in Figure [9] are essentially the same as the previous ones obtained by the 
full approach. All the main features discussed in Section 13.21 are retained. 



5. Conclusions and Discussion 

In conclusions, we develop an approach to extrapolate the coronal magnetic field 
from vector magnetograms based on the Principle of Minimum Dissipation Rate (MDR). 
Analogous to, but yet opposed to the principle of minimum energy, the MDR yields a 
generally non-force free magnetic field through a different variational approach. The 
full MDR-based approach requires two layers of vector magnetograms, while a practical 
approach, representing a class of special solutions to the MDR system, requires only one 
layer, which is more amiable for practical applications to currently available data. A test 
case study using numerical simulation data shows that both approaches recover the solution 
to a good degree of accuracy, as measured by a set of quantitative measures. The errors in 
energy estimate are both within a few percents. Moreover, the non-force free features of the 
solution, mainly the strong perpendicular current density concentrations, are well retained 
in the extrapolation results as well. 

We hereby provide an alternative approach that is fast, and easy to implement, and 
allows one to extrapolate coronal magnetic field in a more general non-force free state 
with manageable effort in a solo work. For instance, all the calculations reported herein 
(128x128x63 grid) were performed on a single-processor 2.8GHz PC in IDL within a 
reasonable time frame. It potentially can be made much faster. It is apparent that the 
Eq. (jll) includes solutions to LFFF (V x B = aB), for —ai = = a. It is not explicit 
whether this is also the case for the nonlinear force-free field (NLFFF) when the parameter 
a i s allowed to vary, alth ough several numerical experiments indicated that this might be 



so (IHu fc Dasgupta 



20071 ) sometimes. 
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We are fully aware of the limitations of the method, due to the fact that the currently 
employed LFFF solver originates from an ill-posed problem. It is adversely affected by 
the increasing resolution of vector magetograms since oimax is inversely proportional to 



N. T 



lis severely limits the search do main of that will yield unique physical solutions 



[e.g., 



Alissandrakis 



extrapolations (e.g., 



981 



Garv 



1989 



Song et a 



Abramenko fc Yurchishin 



Recent progresses in the algorithm for NLFFF 



20061 ). in addition to the existing LFFF algorithm (e.g.. 



19961 ). that modified the problem into a well-posed one, look 



promising and may be adapted. 

It is worth noting that our approach is wholly applicable to extrapolating the plasma 
flow field, by simply replacing magnetic field B with flow vorticity cJ. However the boundary 
conditions are probably harder to obtain since it is the vorticity that is involved. A good 
start is to utilize the flow field data contained in the same 3D MHD simulation test case. 
It is not straightforward to obtain the plasma state since the flow field has to be resolved 
as w ell, and to be consistent, the full set of MHD equations has to be considered (see. 



e.g. 



Montgomery fc Phillips 



19881 ). One may also foresee that extension of this approach 
to the whole solar surface data is of great significance. Further theoretical investigation is 
underway. 



We are grateful to Drs. D. Shaikh, M.S. Janaki, and R. Bhattacharyya for useful 
discussions. HQ and BD acknowledge NASA LWS grant NNX07A073G for partial support. 
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Fig. 1. — Contours of En over the (cti, a^) parameter space for the full MDR-based approach. 
The unlabeled contour lines are of increment 0.01. The dotted lines mark the limit where 
|«i,3| = c^max = '2n / N . The plus sign denotes the location where En is minimum. Due to 
the interchangeability of ai and a^, the plot is symmetric about the dashed line, ai — a^. 
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Fig. 2. — The transverse magnetic field vectors at z = 0. Left: The exact distribution from 
numerical simulation data. Middle: The derived optimal result from the full MDR-based 
approach. Right: The result from a corresponding potential field extrapolation. 
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Fig. 3. — 3D field line plots of the exact solution (left panel) and the corresponding extrap- 
olated solution (right panel) from the full MDR-based approach. All the field lines originate 
from the same set of root points. The bottom image represents the normal magnetic field 
component distribution on the bottom boundary, with the gray scales ranging from strongly 
negative (black) to positive (white). 
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Fig. 4. — Contour plots of field-line-integrated current densities from the exact solution. 
Left The component parallel to the magnetic field. Right: The component perpendicular to 
the magnetic field. The gray colors show the levels of the contour lines, as indicated by the 
colorbars. 



- 19 - 




Fig. 5. — The corresponding field-line- integrated current density distributions obtained by 
the full MDR-based approach. Format is the same as Figure HI So are the gray scales. 




Fig. 6. — The £"„ distribution for the practical MDR-based approach. Format is the same 
as Figure [H 
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Fig. 7. — The corresponding transverse magnetic field vectors at 2; = 0, obtained by the 
practical approach. 




Fig. 8. — The corresponding 3D field line plot, obtained by the practical approach. Format 
is the same as Figure O 




Fig. 9. — The field-line-integrated current densities obtained by the practical approach. 
Format is the same as Figure HI 



Table 1: Figures of merit (see 



Schrijver et al 



(120061 ) and text for definitions). 



Case BP 




Ccs 






e 




Our result 1*^ 


0.97 


0.91 


0.74 


0.62 


0.99 


1.18 


Our result 2^ 


0.91 


0.76 


0.53 


0.27 


1.12 


1.34 


Our result S'^ 


0.97 


0.93 


0.76 


0.64 


0.95 


1.14 


Exact 


1 


1 


1 


1 


1 


1.20 


Potential 


0.95 


0.92 


0.70 


0.63 


0.84 


1.00 



"Results obtained by the full MDR approach (see Section [3]) , where two bottom layers of vector magnetic 
field data are utilized to derive the bottom boundary conditions. 

^Results obtained by the reduced bottom boundary data (see Section [221) , where only the vector magnetic 
field on the bottom layer, and the normal magnetic field component on the layer immediately above are 
utilized. 

^Results obtained by the practical approach (see Section 2]), where only the magnetic field vectors on the 
bottom boundary are required. 



